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In a previous paper, the convergence of a new effective field theory and 
density functional theory (EFT/DFT) approach to the description of the 
nuclear many-body system was studied. The most sophisticated parameter 
set (here Gl) determined by Furnstahl, Serot and Tang from a fit along the 
valley of stability was found to provide quantitative predictions for total 
binding energies and single-particle and single-hole binding energies, spins, 
and parities for selected doubly-magic nuclei far from stability. Binding 
energies of all the even Sn isotopes (Z=50) were also well described. Here 
the calculations are extended to all isotones with magic numbers N=28, 50, 
82, 126, and isotopes with Z=28, 50, 82 with similar predictive results. 

PACS numbers: 21.10.-k; 21.10.Dr; 21.10.Pc 

1. Introduction. 

In a previous paper the convergence of a new approach to the de- 
scription of the nuclear many-body system, based on principles of effective 
field theory (EFT) and density functional theory (DFT) 2J, was studied 
by applying it to selected nuclei far from stability. The results of the cal- 
culations showed an overall agreement with experiment below 1% for the 
total binding energy of the doubly magic nuclei 1 |oSn 8 2, 1 5oSn 50 , 28^20 an d 
2§Ni5o as well as for the entire range of known even-even Sn isotopes. In 
addition, calculations for the chemical potential of nuclei differing by one 
particle or one hole from these doubly-magic ones, produced results with an 
overall agreement with experimental values below the 10% level, and also 
correctly predicted the spin and parity of these nuclei. It was also shown 
in that the calculations converged to the experimental values as the level 
of approximation of the effective lagrangian was increased. From the various 
levels of approximation used, it was found that the ones with the highest 
levels reproduced better the properties of nuclei far from stability. In addi- 
tion to this, it was also observed that when this approach was tested to the 
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whole range of even-even Sn isotopes the results obtained agreed equally 
well with the experimental values. 

Based on the agreement of these results, a calculation of /3-transition 
rates was performed for nuclei close to 1 fQSng2 0- Here the objective was 
to test, through a semi-leptonic process, the accuracy of the single-particle 
Dirac wave functions obtained in describing a particle or a hole out- 
side the filled core. The calculations where done using a current-current 
interaction hamiltonian in which the lepton current is treated in first-order 
perturbation theory while the nuclear currents are expanded in multipoles 
and treated to all orders. The nuclear electroweak currents used in these 
calculations are those derived directly from the same effective lagrangian 
describing the nuclear many-body system. In this way the calculations 
are done based on a self-consistent approach. Two different cases where 
explored in particle-particle transitions fgSns3 — > 1 fiSb82 + e~ + V e 
and hole-hole transitions fglnsa— > foSngi + e~ +V e , The results obtained 
indicated that the particle-particle transitions proceeding between ground 
states were well reproduced (within the 5% level) while transitions to or 
from excited states deviate systematically from the experimental values. 
This behavior was expected based on the fact that ground-state observables 
(in this case the ground-state density and therefore presumably the wave 
functions) are best reproduced on the basis of DFT. Results obtained in the 
hole-hole case, in which all transitions start or end on an excited state, are 
less reliable in reproducing the /3-transition rates. In any case, more studies 
are underway to understand this behavior. 

Both papers showed that this EFT/DFT approach can be used success- 
fully to describe some ground-state properties of the selected nuclei far from 
stability; this extends the applicability of the approach beyond the valley 
of stability. To further explore the range of validity of this approach, the 
calculations done in PP are extended here to all isotones with magic num- 
bers N= 28, 50, 82, 126 and isotopes with Z= 28, 82, 126. Extensions of 
the calculations made in [3] are considered in another paper. 

The calculations done in pP, which are extended here, are based on 
an energy functional of the ground-state density E[p g-Sm ]. By minimiza- 
tion of this energy functional with respect to the density one obtains the 
ground-state energy and density of the nuclear system. The existence of a 
unique functional of the density is guaranteed by the Hohenberg-Kohn the- 
orem |IJ 13] . This minimization process can be computed using a procedure 
developed by Kohn and Sham in which the ground-state density of a inter- 
acting fermion system is computed using the single-particle wave functions 
obtained from an equivalent systems of non-interacting fermions subject to 
an appropriate external potential, known as the Kohn-Sham potential j3J[S|- 
The basis for this procedure is contained in the Hohenberg-Kohn theorem 



eft-dftextension printed on February 8, 2008 



3 



which can be stated by saying that there is a one-to-one correspondence be- 
tween the ground-state density of a system of interacting fermions subject 
to an external potential and the external potential 1 . 

Although the formulation based on DFT is straightforward, the exact 
energy functional is unknown, particularly in the nuclear system where the 
exact form of the nuclear interaction is still a subject of intense research. 
To overcome this, an energy functional is constructed based on the princi- 
ples of effective field theory (EFT) The starting point is an effective 
lagrangian density that satisfies the symmetries imposed by QCD, which is 
considered the underlying theory of the nuclear dynamics. This lagrangian 
is built using the lowest lying hadronic degrees of freedom ( as well as the 
electromagnetic field). The result is the most general lagrangian contain- 
ing (in principle an infinite number of) interaction terms constructed using 
these fields and subject to the imposed symmetries. The various interac- 
tion terms in this lagrangian are ordered in powers of the fields and their 
derivatives. This expansion is equivalent to an expansion in the ratios of the 
fields to the nucleon mass and the Fermi momentum to the nucleon mass. 
These expansion parameters are small: ~l/3 and ~l/4 respectively. Naive 
Dimensional Analysis (NDA), which is a prescription by which the various 
interaction terms are suppressed by powers of a relevant mass scale A (the 
nucleon mass M in this case), indicates that the remaining coupling con- 
stants should be of order unity, a characteristic called "naturalness." This 
natural size of the constants can only be verified after fitting the theory to 
experimental data. This procedure was done in |2j where it is shown that 
the remaining constants are indeed "natural." There the constants were fit 
to ground-state properties of selected spherical nuclei along the valley of 
stability. The effective lagrangian, as well as the corresponding equations 
of motion obtained from it, are shown in the next section. 

As has been mentioned before, this paper explores the extent of appli- 
cability of this EFT/DFT approach. To understand the objectives of this 
work, one has to keep in mind that the coupling constants of the effective 
lagrangian have been fit to properties of selected stable nuclei lying in the 
/3-stability region. Therefore the results of pQ, and the ones presented here, 
test how well this theory reproduces ground-state properties outside the 
stability region. As indicated in the relevant observables which can be 
reproduced, based on DFT, are the ground-state binding energy and chem- 
ical potential of the nuclear system, so this paper concentrates on these 
quantities. The calculations focus on nuclei with semi-magic number of 
neutrons, i.e. N = 28, 50, 82 and 126 or protons, i.e. Z = 28, 50 and 82. 
In general it is seen that the predictions of the total ground-state binding 

1 In relativistic theories the energy functional is written in terms of more than one 
density: e.g. scalar densities 'tpip, vector densities V»7 M ?/>, etc. 
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energy, for both isotones and isotopes, lie below the 1% level as was found 
in pP, and that as the size of the nucleus increases, the error is reduced. 
The agreement between the results of the chemical potential with experi- 
mental data is also consistent with the conclusions reached in pQ. In the 
case of the prediction of the ground-state J n values of odd nuclei neighbor- 
ing even-even ones, the agreement with experimental values varies between 
isotones and isotopes. When the number of neutrons gets very large in the 
regime N>Z, the level density gets very high and the detailed ordering of 
the levels between mayor shells depends sensitively on the details of the var- 
ious densities. This small energy shifts (~ tens keV) make the predictions 
of the ground-state quantum numbers less accurate, though the chemical 
potentials are still relatively well-described. 

All nuclei selected for this study are essentially spherical jS], therefore 
the approximation of spherically symmetric nuclei, made here, is assumed 
to be an accurate one. The calculations are done using the most convergent 
parameter set obtained in [2], called Gl, which corresponds to the most 
sophisticated effective lagrangian. 

To summarize, this paper extends the calculations made in [Q to known 
semi-magic isotones and isotopes using the Gl parameter set. In the follow- 
ing section the formulation of the Kohn-Sham equations and potentials is 
described, and this is followed by a section on results. The paper ends with 
the conclusions and proposes new directions of investigation to further test 
this approach. 

2. Formalism. 

In this section an overview of the formalism of EFT/DFT is presented. 
A more thorough presentation can be found in 

The energy functional that describes the dynamics of the nuclear sys- 
tem is obtained from an effective lagrangian. This is constructed using the 
relevant lowest lying hadronic degrees of freedom which are the pion 
and the rho fields (isovectors) and the nucleon field (isospinor). In addition 
to these, a scalar field is included to reproduced the known short-range at- 
traction of the nuclear force, and a vector field (associated to the u meson) 
that reproduces the mid-range repulsion of the nuclear interaction is also 
included. 

The pion field is introduced as the phase of a chiral rotation of the 
identity matrix in isospin space. A field £(x) is defined in the following way: 

f(a?) = exp(ivr(x)// 7r ). 

where the pion field ir(x) is defined here by ir(x) = \r ■ jr. 
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506.7 


Co 


3.5249 


gi 


97.39 


v P 


-0.2722 


g v 


147.09 


a*i 


1.8549 




77.033 




1.788 




0.0706 


fv 


0.4316 


m 


-0.96161 




4.1572 


^3 


2.2067 


ft 


0.02844 


K4 


-10.090 


Pv 


-0.24992 



Nucleons are included as a two-component isospinor field N(x). The 
upper and lower components correspond to the proton and neutron fields, 
respectively. 

N = ( P ±\ ) (D 



n(x) 

The jOu-field has been added to the lagrangian to include contributions 
coming from the asymmetry between the number of protons and neutrons. 
In addition to the above fields, an electromagnetic field is included and 
the electromagnetic structure of the nucleons is also taken into account. 
Their anomalous magnetic moments, labeled as X p and A n for proton and 
neutron respectively, then enter. 

The effective lagrangian contains, in principle, an infinite number of 
terms. In order to make the theory predictive, these terms have to be ar- 
ranged according to some ordering scheme that takes into account all equally 
contributing terms. Here the various interaction terms in the effective la- 
grangian are arranged following the prescriptions of NDA T . As mentioned 
before, all calculations performed here will be done using the parameter set 
Gl which corresponds to the most sophisticated lagrangian density of 0, 
and is reproduced in Tabled The pion fields are omitted from the follow- 
ing lagrangian since they do not develop a mean field in the class of nuclei 
considered here. 

The nucleonic part of the lagrangian density is 2 : 

C N (x) = N{i^D^-M + g s <P)N- f ^Np^N- f -^Np^N 
-^FpuNXa^N - ^NlpWs + f3 v r 3 )Nd u F^ (2) 



The metric and conventions used are those of [5J|Sj, and differ from |10|. 
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where = + ig p p,j, + igyV^ + §e^ M (l + r 3 ). 
The mesonic part is the following: 



^Co9 2 v(V,Vn 2 + (1 + V P 9 -^)m 2 p tr(p^ 

Ji2 



In the above lagrangian, (ft is the scalar field, the vector field, p^ the 
isovector field rho defined by p^ = \r ' • p^ and A p is the electromagnetic 
field. 

From the above effective lagrangian density one obtains the Kohn-Sham 
equations and the auxiliary equations for the potentials. To do so, the 
various mesonic and electromagnetic fields as considered as local classi- 
cal fields. Then the necessary equations are obtained by calculating the 
Euler-Lagrange equations. Since the nuclei considered here are (essentially) 
spherical jS], the various vector fields: V p , p M and A^, develop only their 
time-components Vq, po and Aq respectively. In addition, the isovector 
po(^)-field develops only its zero-charged component denoted here as bo . 

The Dirac hamiltonian for the nucleon fields as derived from the effective 
lagrangian takes the form: 

h(x) = -ia-V + W(x) + ~T 3 R(x)+P(M-$(x)) + ~(l+T 3 )A(x) 



- ' (fp^nVR + f v VW) + ^p(Ps + P v r 3 )V 2 A 

-—Xpa-VA (4) 
2M v ' 

Here A = ^A p (l + t 3 ) + ^A n (l — t 3 ) and the numerical values used for the 
anomalous magnetic moments are A p = 1.793, A n = —1.913. 

The mean meson and electromagnetic fields are denoted by W = g v Vo, 
$ = 5s0o, R = g P bo, and A = eAg respectively. The quantities f p , f v , (3 s ,f3 v 
are parameters fit to experiment. 

The equations describing the various meson and electromagnetic fields 
are the following: 

^9^ 9^ 9 / n m -? JU 9/ K 3 K 4 $ , p 2 , $ ,m^ TI7 , 

-V 2 $ + m 2 <I> = qips(x) -3> 2 (— + — ) + — (V1 + V2 — )^-W 2 
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(5) 

Here p s is the baryon Lorenz scalar density and g s , m s , K3, K4, rji, 772, 9v> 
g p , ati, «2 are again parameters fit to experiment. Furthermore, 

-V^W + mlW = fi[p B ( x ) + ^V-(j>l(x)r)]-( m + ^)^mlW 

(6) 

Here ps is the baryon density, p\ is the baryon tensor density, p c ^ g is the 
charge density and f v , Co are parameters. In addition, 



-V 2 R + m 2 p R = ^;[P3W + ^V.(^(x)f)]-,^ m ;R-^ fttg (x) 

(7) 



Here p^ and p% are the isovector densities, f p , rj p , g p are parameters and 
<7 7 =5.01, is the coupling of the photon to the w-meson. Finally, 

- V 2 A = eVhg(x) (8) 

where p c h g is the charge density. All parameters used in these equations are 
tabulated in Tabled 

From the above Dirac hamiltonian one can write down the energy eigen- 
value equation which will give the single-particle wave functions that are 
used to construct the ground-state densities. The various mesons and elec- 
tromagnetic fields play the role of the Kohn-Sham potentials. The whole 
set of equations - and Q, below, form the Kohn-Sham equations. 

hip a (x) = E a ip a (x) (9) 

To obtain the single-particle wave functions from © , they are expressed 
in terms of Dirac spherical wave functions: 



M ' ) -[-W$Z.) < ' (10) 



where the <3? Km is a spin spherical harmonic. £t is a two-component spinor 
and the index t is equal to 1/2 for protons and -1/2 for neutrons. By 
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inserting this form of the solution into Eq.JHJ) and using Eq.Q a set of 
two coupled first-order differential equations for the G and F functions is 
obtained: 



ii + r) Ga{r) ~ [Ea ~ Ul{r) + U ^ r)]F ^ r) - UsGa( - r) = 

-) F a (r) + [E a - C/x(r) - U 2 {r)]G a {r) + U 3 F a {r) = (12) 
r J 



d 

s dr 

where the single-particle potentials are given by 



U^r) = W(r) + t a R(r) + (t a + i)A(r) 

+^{f3s + 2t a p v )V 2 A{r) (13) 

U 2 (r) = M-$(r) (14) 
1 

2M 

+ A'(r)[(Xp + A n )/2 + t a (A p - A n )]} (15) 



U 3 (r) = —{f vW '(r) + t a f p R'(r) 



Here the prime indicates a radial derivative; e.g. W (r) = dW(r)/dr. Once 
the single-particle wave functions are calculated, the various densities that 
appear on the r.h.s. of the meson equations can be obtained. They are 
defined as follows: 

OCC r, ■ . i 

=£^(^M-F a 2 (r)) (16) 



a 

OCC 



PbV) = E^(^W + ^(r)) (17) 
^(*) = E^^2G a (r)F a (r) (18) 



Q 
OCC 



a 

occ „ . 

= E^I^( 2 *«) 2G! »(0^(r) (20) 

where the sum goes over the occupied orbitals. In addition to these, the 
charge density is defined as the sum of two pieces: a direct nucleon charge 
density p<i(x) and the vector meson contribution p m (x) 

Pchg = Pd(x) + p m (x) (21) 
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The direct part is given by 

Pd(x) = Pp (x) + ^Lv-(pZ(x)r) + -L^[/3 s V 2 p B (x) + (3 v V 2 p 3 (x)] (22) 

and the vector meson contribution arising from the coupling of the neutral 
vector mesons to the photon takes the form 

p m (x) = —V 2 R + —^—V 2 W (23) 
Here the point proton density p p and nucleon tensor density p^ are given 

by 

P P = ^(PB + Ps) (24) 

occ 

p^ = Y,^(x)iX(3a-rMx) (25) 

All the above equations were derived from the effective lagrangian given 
in Eqs. © and ©, and are the same as those used in ^ and reproduce 
those shown in |2j . To solve them an independently developed program was 
written to solve these local, coupled, non-linear equations 0-© and Eqs. 

and (|12|) . where the appropriate densities are defined in Eqs. (16) - 
(20). This program was used to obtain the results in and reproduces the 
results shown in 0. 

The meson equations are solved iteratively using a Greens function 
method. Because NDA guarantees that each additional term on the rhs 
of Eqs. © - © is smaller than the previous one, convergence is both 
expected and obtained. 

The whole system of Eqs. (fT2*|) and (J5J) - © is solved self-consistently 
until a global convergence is reached. 

3. Results. 

The procedure described above to solve the Kohn-Sham equations to 
find the eigenvalues and wave functions is performed using the parameter 
set Gl given in P3H|, an d which is reproduced in Table \T\£or completeness. 
These equations are solved and used to calculate the ground-state total 
binding energy of the selected nuclei 3 . Before starting the discussion of 
the results, it is important to keep in mind that the effective lagrangian on 
which these calculations are based, was fitted to ground-state properties of 



3 The total binding energy of the system follows directly from the effective lagrangian 
and hamiltonian, see |5] 
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Fig. 1. Calculated and experimental values of the total binding energy for even 
N=28 isotones. Energies are positive and measured in MeV. 



selected spherical nuclei in the stability region [2j and that there has been 
no refitting of these parameters. 

The first set of results shown correspond to the total binding energy of 
even-even isotones with total neutron number, N, equal to 28, 50, 82 and 
126. These are shown in Fig. Fig. As can noted from these figures, 
there is a very good agreement between the calculated and experimental 
values, which are taken from HJ. In Fig. |SJ the percentage deviation, or 
error, of the results is shown. The overall agreement can be seen to improve 
as one goes to larger nuclei and, except for the N=28 isotones, the percentage 
error is in general less than 0.5%. 

The following set of figures, Fig. 0- Fig. 13 show the predictions for the 
ground-state energies of nuclei neighboring even-even isotones of given N. 
The energy differences are taken between the odd nucleus and the closest 
even nucleus. These differences in energy, i.e. the magnitude of the chem- 
ical potential, are calculated using the binding energies obtained in this 
EFT/DFT approach and are compared with available experimental data. 
The even-even isotones shown in these plots, and which are located at the 
zero point of energy, correspond to nuclei with (calculated) closed j-shells. 
A filled j-shell has an even number of particles, and J 71 ^ equal to 0+, here 
only particles or holes relative to filled j-shells contribute to the total an- 
gular momentum and parity of the nucleus. This is precisely the argument 
given in the shell model ^U]. Since these calculations are effectively mean- 
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Fig. 2. Calculated and experimental values of the total binding energy for even 
N=50 isotones. Energies are positive and measured in MeV. 
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Fig. 3. Calculated and experimental values of the total binding energy for even 
N=82 isotones. Energies are positive and measured in MeV. 



field calculations, there is no explicit inclusion of pairing d_. Thus states 
corresponding to partially filled j-shells are not included in these plots. 

From these plots one can appreciate the excellent agreement between 
the calculated and experimental values of this energy difference, or chemical 
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Fig. 4. Calculated and experimental values of the total binding energy for even 
N=126 isotones. Energies are positive and measured in MeV. 
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potential. In addition to this, the predicted ground-state quantum numbers 
of the odd nuclei agree with the experimental values. The only disagreement 
in all the isotone cases considered occurs for ^Tbsj, Fig. E Here the 
experimental ground-state corresponds to a (1/2+) state while the predicted 
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even-even nuclei in the middle row. Calculated J T values are included and, unless 
indicated otherwise, are in agreement with those of the experimental ground state. 
Where the experimental value disagrees with the calculated one, or is uncertain, 
it is also included in the figure. 



ground-state is 11/2-. The result is not far from the metastable state (11/2- 
) which lies 50 keV above the true ground-state This is indicated in 
Fig. |H1 where the metastable state has been shifted in order to be seen on 
this energy scale. 

The next set of plots, Fig. EH - Fig. EH show the ground-state total 
binding energy of even-even isotopes with Z= 28, 50 and 82. They are 
equivalent to the plots shown for the isotones. Here again the predicted 
values agree overall with experiment below the 1% level as can be seen in 
Fig. 1131 The agreement between predicted and experimental data improves 
as one moves to larger nuclei. 

The final set of plots, Fig. El - Fig- QUI correspond to predictions of 
the relative binding energy of odd nuclei neighboring even-even isotopes of 
given Z, with respect to the binding energy of these even-even nuclei. As 
in the equivalent plots for isotones, the even-even nuclei correspond to the 
calculated filled j-shells and locate the zero-energy level. 

Figure 1141 shows the results for isotopes of Ni (Z=28). Here one can 
observe that overall, there is a very good agreement with experiment, similar 
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Fig. 7. Same as in Fig. for N=50 isotoncs. 
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to that of the N=28 isotones. In the case of 28M33, there is a disagreement 
in the predicted quantum numbers of the ground-state. The experimental 
value corresponds to a 3/2- state, while the predicted one is a 5/2-. There is 
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Fig. 10. Calculated and experimental values of the total binding energy for even 
Z=28 isotopes. Energies are positive and measured in MeV. 



an excited state with J w value in agreement with the calculated ground-state 
that lies only 70 keV above it In this figure, this separation in energy 
has been exaggerated to show it in the scale of the plot. The disagreement in 
the location of the ground-state shows the sensitivity of the actual ordering 
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Fig. 12. Calculated and experimental values of the total binding energy for even 
Z=82 isotopes. Energies are positive and measured in MeV. 



of states to the form of the central potential. One can understand this by 
looking at the effects of putting additional neutrons into the system. There 
is a major shell closure at N=28. If one more neutron is added to the system 
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even isotopes considered. 



it will go into the next j-shell which is j = 3/2 corresponding to a 2p3/2 level, 
as expected from the naive shell model 4 . This level becomes completely 
filled at N=32, since 2j + 1 = 4, which is correctly predicted as can be seen 
in the plot. Adding one more particle to this nucleus, should start filling the 
next j-shell, the lf5/2 level, giving the ground state the quantum numbers 
5/2- as predicted. Instead, it seems that it is energetically more favorable to 
promote one particle from the 2p3/2 level to the lf5/2 leaving a hole. The 
two particles in the lf5/2 level do not contribute to the angular momentum, 
which is solely determined by this hole in the 2p3/2. Adding more neutrons 
fills the lf5/2 until it is completely full at N=38, since 2j + 1 = 6, which 
again is correctly predicted. After this, more neutrons just start filling the 
next level, 2pl/2, following the ordering given by the naive shell model. 
Therefore the predictions follow very closely the way the neutron levels are 
being filled and the fact that the disagreement at N=33 is small only shows 
how sensitive the ordering is to the form of the central potential. 

The exact ordering of energy levels has a larger effect when the level 
density is high and the levels are very close together. This is precisely what 
happens in the case of isotopes with large values of Z. For any given Z the 
density level of neutrons close to the Fermi surface is larger since in general 
N>Z making the energy levels very close to each other. This can be appreci- 




4 The level ordering is sensitive to the actual form of the central potential, and here 
we use the one presented in |1U| . chapter 6. 
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ated in Fig. ^Jwhere, if there is a disagreement between the calculated and 
experimental ground-state value of J 71 ", the lowest-lying excited state with 
the calculated quantum numbers has also been plotted ^3]. Note in partic- 
ular the small (~ tens of keV) energy differences between these states and 
each case and the fact that there is always a partner with quantum numbers 
that agree with the ones predicted for the ground-state. As for similar cases 
before, the energy difference has been exaggerated to be appreciated at this 
scale. 

A similar effect occurs in the case of the isotopes of Pb (Z=82) . Figure lT^l 
shows the filled j-shell states closest to 2 g|Pbi26 which is a mayor shell 
closure. Adding an extra neutron to this nucleus will start filling the next 
j-shell, which according to the naive shell model is a 2g9/2 state. This 
is precisely what is predicted here as can be seen in this figure for the 
particle state at N=126 (i.e. N=127). Removing a neutron from 2 g§Pbi26 
corresponds to depleting the 3pl/2 state, as can be seen in the figure, and 
again it is also correctly predicted. Since 2 8fPki26 formed part of the nuclei 
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used to fit the parameters of the theory, it is expected that there will be an 
accurate description of the neighboring nuclei, and in fact there is. 

Removing one additional neutron completely depletes the j = 1/2 shell 
bringing the total number of neutrons to N=124, which is also observed 
occurring in this plot. According to ^U] (see footnote above) the next hole 
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level would be a lil3/2 level, i.e. 13/2+, but the experimental evidence 
indicates that this level is in fact a 5/2-. The reordering of the levels, 
with respect to the naive shell model, can be understood as either a large 
splitting between the 3pl/2 and 3p3/2 levels, which brings the 3pl/2 level 
above the 1H3/2, or a considerable lowering of the 1H3/2 which brings the 
1H3/2 below the 3pl/2. In either case removing one more neutron from 
this nucleus will start depleting the next j-shell which could be a 3p3/2 or 
a 2f5/2. In Fig. ^] it is shown that this level is a 5/2- in disagreement with 
the calculated one, 3/2-. The closest excited state having this predicted 
quantum numbers lies only 262 keV above the ground-state. This state is 
also plotted [Tl|- The energy difference has been exaggerated in the figure, 
as in the previous plots. The reordering between the 2f5/2 and 3p3/2, shift 
the position of the real j-shell closures with respect to the calculated ones, 
as happened for the Z=50 isotopes shown in Fig. 1151 therefore they have 
not been included in this plot. Despite this disagreement, the magnitude of 
the chemical potentials of these nuclei are relatively well described. 

In summary, these results show that the binding energy of all even-even 
nuclei with semi-magic numbers of neutrons and protons, N = 28, 50, 82, 
126 and Z = 28, 50, 82 are reproduced to 1% or better in this EFT/DFT 
approach with parameters Gl fit along the valley of stability. With all the 
isotones, the chemical potential of nuclei with a single proton particle or 
hole outside of filled j-shells is well reproduced and the predicted value of 
J 71 ' for the ground-state of these nuclei agrees with experiment. 

When the number of neutrons gets very large in the regime where N 
> Z, the level density is very high and the detailed ordering of the levels 
between major shells depends sensitively on the details of the various densi- 
ties. Thus, while the chemical potential (binding energies) are still relatively 
well-described, the theory begins to lose its predictive power for the spin 
and parity of the ground states of single neutron particle and hole nuclei 
neighboring those with filled j-shells. 

4. Conclusions. 

This work investigates the applicability of the EFT/DFT approach to 
reproduce total binding energies and chemical potentials of semi-magic nu- 
clei (N = 28, 50, 82, 126 and Z= 28, 50, 82) far from the valley of stability. 
It extends the work done in d], where the convergence of this EFT/DFT 
approach was studied for doubly- magic nuclei. All calculations were per- 
formed using an independently developed program to solve the Kohn-Sham 
equations derived from the most sophisticated effective lagrangian given 
in [2] with parameter set Gl. 

From the figures given in the previous section, several conclusions are 
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drawn: 

1. The EFT/DFT approach can reproduce the total binding energies of 
even-even isotones and isotopes with semi-magic numbers of neutrons 
and protons, respectively. The agreement is in general better than 1 % 
where the best results are obtained for larger nuclei. This conclusion 
is in agreement with the previous results obtained in [Q. 

2. The predictions for the magnitude of the chemical potentials for odd 
nuclei corresponding to particle or holes neighboring filled j-shells 
agree overall with the experimental evidence. This agreement varies 
for different nuclei and on average lies between 10% and 20%. 

3. Regarding the predictions of the ground-state quantum numbers of 
odd nuclei neighboring even-even nuclei with filled j-shells, the results 
indicate that in the case of isotones these are accurately reproduced; 
for heavy isotopes with N > Z, although the predictions of chemical 
potentials are in essential agreement with experiment, the EFT/DFT 
approach begins to lose its predictive power. 

Overall, the use of this effective field theory approach for QCD and 
the interpretation of the resulting relativistic Hartree equations through 
Kohn-Sham potentials and density functional theory [2] provides an excel- 
lent description of the ground-state properties of semi-magic nuclei when 
extrapolated away from the valley of stability. 
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